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Abstract 

Positronium (Ps) reaction rates (k) with weak Acceptors (Ac) leading to the formation 
of Ps-Ac Complexes show several intriguing features: non-monotonic temperature depen- 
dence of K (departing from the usual Arrhenius paradigm), considerable variability of k 
with respect to different solvents, and anomalies in response to external pressure at ambi- 
ent temperature (large changes of k in some media and hardly any in others). We explain 
all these phenomena, introducing the novel concept of a critical surface tension, which 
unifies observations in diverse non-polar solvents at different temperatures and pressures. 

The Positronium (Ps) atom, as the 'lightest isotope of hydrogen', enters into various chemical 
reactions, amongst which we shall concentrate here on those with diamagnetic organic com- 
pounds (equipped with a suitable low-lying molecular orbital) such that a Positronium- Acceptor 
(Ps-Ac) Complex can form in the presence of some solvent (S). 

The observed rate constant (k) for a given Acceptor depends on the solvent and for a given 
solvent varies in a rather remarkable manner with temperature (T) , namely: 

• K increases with T at low temperatures, 

• K, reaches a maximum at T = Tq (the turnover temperature), 

• K decreases with increasing T above Tq. 

This is shown in Fig.l for the case of the weak Acceptor nitrobenzene in different solvents. 
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Fig.l Observed temperature dependence of rate constants for Ps reaction with nitrobenzene in 
different solvents. (Taken from W.D. Madia et.ai, J.Am.Chem.Soc. 97(1975)5041)- The lines 
correspond to slopes in the enthalpy dominated and diffusion controlled regimes as predicted by 
our model(with appropriate over-all normalization for heptane). 

This unusual phenomenon [1], first observed by Goldanskii and his group [2] at Moscow, 
and subsequently by Hall, Madia and Ache [3], is in sharp contrast with the normal trend in 



activated chemical processes where In k versus 1/T (the inverse temperature) yields a straight 
line with negative slope, as the rate is proportional to the Maxwell-Boltzmann probability 
factor exp{—E*/kBT) where ks is the Boltzmann constant and E* is the activation energy or 
barrier height to be overcome by thermal agitation. Possible prefactors are generally taken to 
be mildly dependent on temperature. 

Another dramatic observation [4,5] pertains to the variation of the rate constant for Ps-Ac 
Complex formation with external pressure at ambient temperature, which shows surprisingly 
strong solvent dependence. Thus, for instance, with nitrobenzene as the Acceptor and hexane 
as the solvent the rate constant k, was found to be enhanced by a factor of almost thirty as 
the pressure was increased to ~ 1000 kg/cm^, while with benzene as the solvent the same rate 
manifested only small variation over the same pressure range. This feature is depicted in Fig. 2. 
Kobayashi [4,5] conjectured that this could indicate a difference in reaction mechanism in cases 
which do and those which do not show marked pressure dependence, and that a correlation 
could be sought with the corresponding solvent dependence of the turnover temperature. 
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Fig.2 Observed dependence on external pressure for Ps reaction with nitrobenzene in hexane 
and in benzene. The predicted dependence in the case of hexane is shown by the solid line 
(appropriate normalization has been done). 



The objective of this letter is to put forward a heuristic model which captures the essence 
of all these puzzling aspects of Positronium reactions outlined above, through a remarkably 
simple and reasonably accurate semi-quantitative description. As we shall see, the widely used 
bubble model for Ps-annihilation in liquids provides a suitable setting for the purpose at hand. 
This was proposed by Ferrel [6] to account for the observed largeness of the ort ho- Positronium 
(o-Ps) pick-off lifetime, whereby the positron in o-Ps seeks out an electron in the surrounding 
medium with opposite spin (not its partner in o-Ps which is in a spin triplet state) to decay into 
two gammas. He suggested that due to the repulsive interaction between o-Ps and the atoms 
of the surrounding liquid, the Positronium pushes away the molecules of the liquid and gets 
self-trapped in a cavity or bubble. The decreased probability for the positron to find an electron 
(of the surrounding medium) in its immediate vicinity leads to the lengthening of the pick-off 
lifetime. For simplicity he took the trapping potential to be an infinitely repulsive spherical 
well of radius R. The resulting quantal zero-point energy of Ps with mass 2m ( m being the 
electron mass) is Eq = ^^^2 ■ This exerts an outward 'pressure' viz. — to be balanced by the 
contractile forces of compression due to the surface tension (a) of the fluid, thereby minimizing 



the total energy Etot = Eq + AttR'^o- viz. 



a+Syro" = resulting in a bubble of radius 



Rq = {^f/^ = (12.445)/ai/4A , where a is measured in dynes/cm ( .624x 10 eV/A^). 



This model was further developed by Tao [7] and by Eldrup et al. [8] employing again the 
infinite spherical well but describing the picked-off electrons as forming a thin layer uniformly 
coating the inner surface of the bubble. Though such a description has been criticised on 
account of the infinite repulsion [9,10] or the unrealistic nature of the sharpness of the bubble 
edge [11,12], nevertheless as a first approximation the model does rather well in explaining the 
observables in the case of pure solvents (lifetime and angular correlation of decay gammas). 

In this framework the formation of the Ps-Ac Complex within the cavity results, through 
the restriction on the motion of the Positronium due to its binding to the Acceptor, in a release 
of the outward 'uncertainty pressure' and accordingly the bubble begins to shrink under the 
infiuence of the surface tension. As the bubble wall approaches the Positronium (bound to the 
Acceptor with some energy E = —B where B is the binding energy), the increasing proximity 
of the solvent molecules leads (because of the augmented repulsive energy) to a reduction of 
the binding, until at some radius Rh of the bubble the Ps-Ac bond is broken viz. 5 = 0. To 
raise the status of this contention [13] to at least a semi- quantitative level let us consider the 
Positronium to be subjected to an attractive potential V{r) due to the Acceptor, while it also 
feels the infinite repulsion due to the solvent located at the bubble wall at a distance R. Taking 
the potential V{r) to be a spherical well [14] of depth Vq and range a viz. — V^)©(a — R) where 
© is the step function, the relevant Schrodinger equation for the centre-of-mass motion of Ps is 



4m dr"^ 

Am dr'^ 



Va 



u — —Bu for r < a (la) 



u — —Bu for r > a (16) 



with ip = u/r. The appropriate boundary conditions are u{r = 0) = and u{r = R) = 0, while 
u and its derivative must be continuous at r = a, yielding the eigenvalue condition 

k tan Ka — —k tanh k{R — a) (2) 



where k = y and k — \/g^ — k'^ with ^^3^ = g"^- The bond breaking radius Rj,, viz. where 
A; = 0, is thus given by 

tan^a 



Rb^al — . (3) 

L ga J 

Since in the leading approximation (adopted by us) the effect of the solvent is represented by 
an infinite repulsion (tantamount to a boundary condition forcing ijj to vanish at R), the radius 
Rb depends only on the characteristics of the Acceptor vis a vis its affinity for the Positronium, 
and is, to that extent, same for all solvents. However, we must also consider the fact that 
the Ps-Ac-bubble system would equilibrate at some radius Rmin at which the total energy is a 
minimum, namely, 

±[-B + 4nR'a]^0. (4) 

This leads to 

1 1 6C^ 

Rmin = (^^(^'^i-^^ _ tanC)(3C - 3tanC + 3Ctan2 ( - 2tan3C) 

where ( = ga. Whether the Ps- Ac-bubble system will be stable or not will depend on whether 
Rmin is greater than or lesser than the bond-breaking radius Rb. 

A big step forward is made possible by our recognition that in leading approximation Rh 
is essentially solvent independent, while solvent dependence enters through Rmin and that 
too only through a. As the surface tension of the solvent is a function of temperature, the 
equilibrium radius of the bubble (with Ps-Ac inside) depends on temperature only through a 
viz. RminiT) = RminW{T)]. As T is decreased, a increases and by virtue of eq.(4) and more 



explicitly eq.(5), Rmin decreases. When Rmin becomes less than Rf, the Ps- Ac-bubble system 
is no longer stable. Clearly there is a critical value of a, say acr, at which Rmin[ccr] — Rb, 
marking the watershed beyond which the Ps-Ac-bubble system destabilises. The value of T at 
which a attains the value dcr depends on the particular solvent, but acr does not, but is, in 
the leading order, a property of the Acceptor under consideration. We argue that the turnover 
temperature Tq is the temperature at which, for the solvent at hand, cr(To) = acr- Indeed at 
sufficiently high temperature where a{T) < (Jcr, a large negative change in enthalpy occurs 
as the reactants [Ps in a bubble of radius Rq — (jf^)^'^^ and the Acceptor] react to form 
the product [Ps-Ac Complex in bubble] with the radius having shrunk to Rmin- This makes 
it a down-hill reaction, since for weak Acceptors the activation energy is expected to play a 
sub-dominant role and also the effect of solvent viscosity is negligible (because the complex is 
protected by the bubble from the buffetting by the solvent molecules) . It is thus basically the 
negative activation volume that is responsible for the anti-Arrhenius behaviour for T > Tq. On 
the other hand for T < Tq when a{T) > acr the Ps-Ac-bubble system is unstable, the Ac is 
squeezed out of the bubble and the formation of the Complex must of necessity take place in the 
milieu of the solvent accompanied by the continual impact of the solvent molecules. As such, 
following Kramers [15] the role of the medium would be said to belong to the Smolochowski 
regime, with the dependence of the rate constant on the viscosity (77) varying as ^ ~ D (where 
D is the diffusion coefficient of the liquid, the last step being a consequence of the Einstein- 
Stokes relation). In view of the smallness of the activation energy for weak Acceptors, we would 
therefore expect that k. ~ exp[—E^/kBT], where is the activation energy associated with 
the process of diffusion. The reaction in this region would be diffusion controlled and would 
exhibit a normal Arrhenius behaviour. This enables us to make the important prediction that 
a(To) — acr and that it is approximately solvent independent. In order to confront this with 
experiment, we have plotted in Fig. 3, a{T) against T — Tq for various solvents (for which data 
is available with nitrobenzene as the Acceptor). 
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Fig.3 Surface tension a{T) for different solvents as a function ofT — TQ exposing the concept 
of the critical surface tension (acr) 

It is indeed highly gratifying to note that while Tq differs widely from solvent to solvent, 
and also the values of the surface tension cr at a given temperature for different solvents have a 
substantial spread, nevertheless o-(To) = for the solvents under consideration lie in a rather 
narrow range, namely, acr = 26 ± 2 dynes/cm. 

Armed with this value of acr we are now in a position to access the Ps-Ac interaction 
parameters which, as we shall discuss below, are unfortunately not available from any other 



source. Being led by the estimates of other authors [16] let us fix the range of the interaction 
to be a = 1.5 A , and determine g by putting, in accordance to our discussions above, a — acr 
and Rmin = Rb [refer eq.(3)]. This immediately yields the value g =1.25 which corresponds 
to the Ps-Ac binding energy (the Acceptor being nitrobenzene) Bq = O.lSeV. Note that Bq 
signifies the basic binding in the absence of the solvent, namely, Bq = Limji^ooB , which implies 
vide eq.(2), the well-known result: —k cot kg = k the eigenvalue condition for a spherical well. 
Unfortunately there is no direct measurement of this binding energy Bq, and the approaches 
based on first principles [17] are beset with huge theoretical uncertainties. However, using our 
methodology we have now completely tied down the model and shall proceed to show how this 
explains all the main observed features, 

In order to find the slope of the In k versus 1/T plot in the higher temperature anti- 
Arrhenius region we recognise that this is nothing but the activation free energy, the major 
contribution to which come from the change in Enthalpy arising from the shrinkage of the 
bubble from its value given by Rq to the size determined by Rmin, viz. AnRfj^^^cT — A-kRqCt 
where, of course the radii are functions of temperature through the surface tension. To this 
must be added the change in the Ps-Ac binding due to the approach of the bubble wall from Rq 
to Rmin- The resulting behaviour (appropriately normalized) is depicted in Fig.l for the case of 
Ps nitrobenzene reaction in heptane, for instance, as a solvent. The slope obtained through our 
model corresponds to ~ 0.15 eV which agrees very well with the experimental value. The part 
of the plot in the Arrhenius regime is also shown by using the value of the activation energy 
for diffusion {E^). Fig.l demonstrates that the main aspect of the temperature dependence of 
the rate has been captured in a satisfactory manner. 

Apart from providing an interpretation of the temperature dependence of the reaction as set 
forth above, this simple model is also able to furnish an explanation of the observed variation 
with the external pressure (P). In view of the small compressibility of liquids one would hardly 
expect any appreciable activation volume when T < Tq as the reaction occurs in the solvent 
itself. However, for T > Tq due to the involvement of the bubble a large AV* (activation 
volume) becomes possible and hence a significant pressure dependence can occur, in view of 
the basic Polyani relationship: 

fdlnK\ AV* 

between the effect of external pressure on the rate of a chemical reaction and the activation 
volume of the reaction (namely, the difference between the volume of the activated complex 
and the volume of the reactants). To obtain quantitative estimates it is to be noted that the 
initial volume viz. Ps in the bubble without the Ac must be found now by minimizing the total 
energy, namely Etot = ^^^2 + 4:71 R^ a + ^R^P with respect to R and thereby obtain Ro{P). 
Here P is measured in kgwt/cm^ (=.613 x 10~® eV/A"^). Similarly the equation determining 
the equilibrium radius [Rmin{P)] with Ps-Ac inside the bubble must be appropriately modified, 
so that in place of eq. (4) , we now have 

^{-B + A7rR'a + fR'p)^0 (7) 

Equipped with Ro{P), Rmin{P) and B[R[P)] the variation of the rate constant with pressure 
can be readily calculated. 

Consider the two typical cases depicted in Fig.2 out of the many taken from the experimental 
studies of Kobayashi [4,5]. All these experiments were performed at ambient temperature 
{Texpt — 19± 1°C). With nitrobenzene as the Acceptor and benzene as the solvent we note that 
there is no significant effect of pressure on the reaction rate. This is quite consistent with the 
model in view of the fact that the turnover temperature in this case [18] is considerably above 
Texpt and hence one is working in the Arrhenius region where the pressure effect is expected 
to be small. On the other hand when the solvent is hexane, the experiment was performed 



in the regime where the Ps- Ac-bubble system is stable (as the turnover temperature here is 
-53°C), and accordingly the reaction rate responds appreciably to external pressure. Indeed 

the enhancement of the rate by a factor of about thirty in our model is even in quantitative 
agreement with the experimental results of Kobayashi. Using the experimental value of the 
surface tension of hexane at Tg^pt (cr ~ 19 dynes/cm), the model indicates the behaviour shown 
in Fig. 2 with appropriate over-all normalization. 

Above a critical pressure the Ps- Ac-bubble system destabihzes and the rate versus pres- 
sure curve levels off. Except for the fact that there is a somewhat precocious onset of the 
Arrhenius regime the general experimental trend is captured very well indeed considering the 
approximations involved. 

We may thus conclude that the simple model we have proposed effectively accounts for 
the temperature, pressure and solvent dependence of Ps-Ac Complex formation reactions. Of 
cardinal importance for our discussion has been the introduction of the notion of the critical 
surface tension (ucr) which enables us to semiquantitatively account for most observations and 
also to make the important qualitative remark: that for any solvent at a temperature such that 
its surface tension is smaller than the reaction rate decreases with increasing temperature 
and is also significantly affected by external pressure, while otherwise it is the other way around. 
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